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Abstract. A search for the Standard Model Higgs boson at the Large Hadron Collider (LHC) running at 
a centre-of-mass energy of 7 TeV is reported, based on a total integrated luminosity of up to 40 pb _1 col- 
lected by the ATLAS detector in 2010. Several Higgs boson decay channels: H -> 77, H -> ZZ M -> 
H ^ ZZ -> Mm, H ^ ZZ ~^ Mqq, H -> WW M -> and # -> WW -> (t is e, /i) are combined 
in a mass range from 110 GeV to 600 GeV. The highest sensitivity is achieved in the mass range between 
160 GeV and 170 GeV, where the expected 95% CL exclusion sensitivity is at Higgs boson production 
cross sections 2.3 times the Standard Model prediction. Upper limits on the cross section for its production 
are determined. Models with a fourth generation of heavy leptons and quarks with Standard Model-like 
couplings to the Higgs boson are also investigated and are excluded at 95% CL for a Higgs boson mass in 
the range from 140 GeV to 185 GeV. 



1 Introduction 

The search for the Standard Model Higgs boson [IH3] is 
one of the key aims of the Large Hadron Collider (LHC) 
at CERN. Prior to the LHC, the best direct information 
is a lower limit of 114.4 GeV, set using the combined re- 
sults of the four LEP experiments _4J, and an excluded 
band of 158 GeV to 173 GeV from the combined Teva- 
tron experiments 5 , 6 J . First results from the ATLAS ex- 
periment are available in various Standard Model Higgs 
boson search channels [7TH1|. There are also results from 
the CMS collaboration [L2] in the H -> fffW -> 
channel which have a sensitivity similar to the equivalent 
search reported here. These results are based on proton- 
proton collision data collected in 2010 at a centre-of-mass 
energy of y/s = 7 TeV. This paper combines the results 
from the different Higgs boson searches to obtain the over- 
all sensitivity to a Standard Model Higgs boson with the 
2010 ATLAS dataset. 

All analyses use the most detailed calculations avail- 
able for the cross sections, as discussed in Section [31 The 
searches in individual Higgs boson decay channels H —> 
77, H -> WW^ and H -> ZZ^> are outlined in Sec- 
tions [H and |6l respectively. The statistical interpreta- 
tion uses the profile-likelihood ratio [T3] as test-statistic. 
Thirty-one Higgs boson masses, in steps of 10 GeV from 



110 GeV to 200 GeV (plus 115 GeV in addition) and 
20 GeV from 200 GeV to 600 GeV, are tested. Exclu- 
sion limits are obtained using the power constrained CL s j 
limit Q3] , as discussed in Section [7] To allow for com- 
parisons with the exclusion limits obtained by other ex- 
periments, the results are also determined using the CL S 
method [T5]. The limits are presented in terms of a/a^M, 
the multiple of the expected Standard Model cross sec- 
tion at the Higgs boson mass considered. Results are also 
presented in terms of the corresponding ratio where the 
cross section in the denominator includes the effects of a 
fourth generation of heavy leptons and quarks with Stan- 
dard Model-like couplings to the Higgs boson. Section [5] 
describes the treatment of the major sources of systematic 
uncertainty in the combined likelihood. The limits for in- 
dividual channels and the combined results are detailed in 
Section |H] and the conclusions are drawn in Section [TU] 



2 The ATLAS Detector 

The ATLAS experiment |16j is a multipurpose particle 
physics apparatus with forward-backward symmetric cylin- 
drical geometry covering \rf\ < 2.5 for tracks and \t) < 
4.5 for jetifl. The inner tracking detector (ID) consists 
of a silicon pixel detector, a silicon microstrip detector 



1 In this paper, the raised index '*' implies a particle off 
mass-shell, £ is always taken to mean either e or /j, and q can 



be any of u, d, s, c or 6. 



2 ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of the de- 
tector and the z-axis coinciding with the axis of the beam pipe. 
The a;-axis points from the IP to the centre of the LHC ring, 
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(SCT), and a transition radiation tracker (TRT). The ID 
is surrounded by a thin superconducting solenoid provid- 
ing a 2 T magnetic field, and by high-granularity liquid- 
argon (LAr) sampling electromagnetic calorimeters. An 
iron-scintillator tile calorimeter provides hadronic cover- 
age in the central rapidity range. The end-cap and for- 
ward regions are instrumented with LAr calorimetry for 
both electromagnetic and hadronic measurements. The 
muon spectrometer (MS) surrounds the calorimeters and 
consists of three large superconducting toroids, each with 
eight coils, a system of precision tracking chambers, and 
detectors for triggering. 

The data used in this analysis were recorded in 2010 
at the LHC at a centre-of-mass energy of 7 TeV. Appli- 
cation of beam, detector, and data-quality requirements 
results in a total integrated luminosity of 35 to 40 pb -1 
depending on the search channel, with an estimated un- 
certainty of ±3.4% [T7]. The events were triggered either 
by a single lepton or a pair of photon candidates with 
transverse momentum (pr) thresholds which were signifi- 
cantly below the analysis offline requirements. The trigger 
introduces very little inefficiency except in one channel, 
H — > WW —> tvqq, where there are efficiency losses in 
the muon channel of about 16%. 

Electron and photon candidates are reconstructed from 
energy clusters recorded in the liquid-argon electromag- 
netic calorimeter. The clusters must have shower profiles 
consistent with those expected from an electromagnetic 
shower. Electron candidates are matched to tracks recon- 
structed in the inner detector, while photon candidates 
require either no track or an identified conversion can- 
didate. Muon candidates are reconstructed by matching 
tracks found in the inner detector with either tracks or hit 
segments in the muon spectrometer. Details of the quality 
criteria required on each of these objects differ amongst 
the analyses discussed here. There are in addition isolation 
criteria which again depend upon the specific backgrounds 
relevant to each analysis. 

Jets are reconstructed from topological clusters [TB] in 
the calorimeter using an anti-k t algorithm |19j with a ra- 
dius parameter R = 0.4. They are calibrated [T51I2T)] from 
the electromagnetic scale to the hadronic energy scale us- 
ing pt and r] dependent correction factors based on Monte 
Carlo simulation and validated on data. They are required 
to have a pt greater than 25 GeV unless otherwise stated. 
B tagging is performed using a secondary vertex algorithm 
based upon the decay length significance. A selection re- 
quirement is set to describe a jet as '6-tagged' which has a 
50% efficiency for true &-jets. The missing transverse en- 
ergy is reconstructed from topological energy clusters in 
the ATLAS calorimeters, with corrections for measured 
muons. 



and the y-axis points upward. Cylindrical coordinates (r, <jj) 
are used in the transverse plane, <j> being the azimuthal angle 
around the beam pipe. The pseudorapidity is denned in terms 
of the polar angle 6 as r] — — lntan((9/2). 



3 Cross sections, decays and simulation tools 

3.1 Search for the Standard Model Higgs boson 

At the LHC, the most important Standard Model Higgs 
boson production processes are the following four: gluon 
fusion (gg —> H), which couples to the Higgs boson via a 
heavy-quark triangular loop; fusion of vector bosons radi- 
ated off quarks (qq — > qqH); associated production with 
a vector boson (qq —> WH/ZH); associated production 
with a top-quark pair (qq/ gg — > tiH). The current calcu- 
lations of the production cross sections have been gathered 
and summarised in Ref. |21| . 

Higher-order corrections have been calculated up to 
next-to-next-to-leading order (NNLO) in QCD for the gluon 
fusion [22-27], vector boson fusion [28] and associated 
WH I ZH production processes [53] , and to next-to-leading 
order (NLO) for the associated production with a ti pair [3H1 
131] . In addition, QCD soft-gluon resummations up to next- 
to-next-to-leading log (NNLL) are available for the gluon 
fusion process [22- The NLO electroweak (EW) correc- 
tions are applied to the gluon fusion [33][34] , vector boson 
fusion 135.36] and the associated WH/ZH production [3"T] 
processes. 

The Higgs boson decay branching ratios used take into 
account the recently calculated higher order QCD and 
EW corrections in each Higgs boson decay mode [2T1I38] . 
The errors in these calculations for the states considered 
here are at most 2% and are neglected. For most four- 
fermion final states the predictions by Prophecy4f [391140] 
are used which include the complete NLO QCD+EW cor- 
rections with all interference and leading two-loop heavy 
Higgs boson corrections to the four-fermion width. The 
H -> ZZ — > liqq and H — > ZZ — > livv analyses use the 
less precise single Z boson decay rates from Ref. [41] . 

The total signal production cross section in pp col- 
lisions at -y/s = 7 TeV, multiplied by the branching ratio 
for the final states considered in this paper, is summarised 
in Fig. [T] as a function of the Higgs boson mass. Sources 
of uncertainties on these cross sections include missing 
higher-order corrections, imprecise knowledge of the par- 
ton distribution functions (PDFs) and the uncertainty on 
the strong force coupling constant, a s . These uncertain- 
ties are treated according to the recommendations given 
in Refs. [2TU52H 55] and are ±(15-20)% for the gluon fusion 
process, ±(3-9)% for the vector boson fusion process and 
±5% for the associated WH/ZH production process. 

3.2 Higgs boson search in fourth generation models 

Models with a fourth generation of heavy leptons and 
quarks with Standard Model-like couplings to the Higgs 
boson enhance its production cross section in gluon fusion 
by a factor of 4 to 10 compared to the predicted rate with 
three generations [3Bl449| . The model considered here [50] 
has very heavy fourth generation fermions, giving a min- 
imum cross section but excluding the possibility that the 
Higgs boson decays to heavy neutrinos. These can weaken 
the exclusion for Higgs boson masses below the W pair 
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Fig. 1. The cross section multiplied by decay branching ratios 
for Standard Model Higgs boson production in pp collisions at 
a 7 TeV centre-of-mass energy as a function of mass [21]. All 
production modes are summed, and only final states considered 
in this paper are shown. Two bands are shown for each curve; 
the inner represents the QCD scale uncertainty and the outer 
also includes the a a and PDF uncertainty. 



threshold [51]. It should be noted that the branching ra- 
tio into photons is suppressed by a factor around 8 in this 
model. 

The Higgs boson production cross section in the gluon 
fusion process and its decay branching ratios have been 
calculated in the fourth generation model at NLO with 
HIGLU [5J] and HDECAY [35]. The NNLO+NNLL QCD 
corrections are applied to the gluon fusion cross sections. 
The QCD corrections for the fourth generational model 
are assumed to be the same as in the Standard Model. The 
full two-loop Standard Model electroweak corrections [33] 
[34] are taken into account. The effect of a fourth genera- 
tion in the Standard Model background processes, which 
includes contributions from loop diagrams, has been ne- 
glected. 



3.3 Monte Carlo simulations 

For the H— >ZZ Monte Carlo samples, the Higgs signal is 
generated using PYTHIA [S3] interfaced to PHOTOS [51] 
for final-state radiation. The H — > WW^ — > Ivlv events 
produced by gluon fusion or vector boson fusion are mod- 
elled using the MC@NLO [551155] and SHERPA [57] Monte 
Carlo generators, respectively. H — > WW — > ivqq is mod- 
elled using PYTHIA for the gluon fusion and HERWIG [58] 
for vector boson fusion. The 77 signal is simulated with 
MC@NLO, HERWIG and PYTHIA for the gluon fusion, 
vector boson fusion and associated production processes 
respectively. 

For background sample generation, the PYTHIA, ALP- 
GEN [5J, MC@NLO, MADGRAPH [60], SHERPA and 
HERWIG packages are employed. 



All Monte Carlo samples are processed through a com- 
plete simulation of the ATLAS detector [61 using the 
GEANT programme [6J]. 



4 Search for H — > 77 

The search for the Higgs boson in the 77 decay mode is de- 
scribed below; further details can be found in Ref. [7, . The 
event selection requires the presence of at least two iden- 
tified photons [53], including converted photons, isolated 
from any other activity in the calorimeter. The leading 
and the sub-leading photons are required to have trans- 
verse momenta above 40 GeV and 25 GeV, respectively. 
The directions of the photons are measured using the po- 
sition determined in the first sampling of the electromag- 
netic calorimeter and that of the reconstructed primary 
vertex. The di-photon invariant mass spectrum is used to 
search for a peak above the background contributions. 

The main background processes in the H — > 77 search 
arise from the production of two isolated prompt photons 
(77) and from fake photons in photon-jet (77) and di-jct 
(jj) events. Fake photons can originate from jets in which 
a leading ir° or r\ meson from the quark or gluon fragmen- 
tation is reconstructed as a single isolated photon. Each of 
these background contributions has been estimated from 
sideband control samples in the data. The backgroud from 
Droll- Yan events, Z/j* — > ee, where the electrons are mis- 
takenly identified as photons, is estimated from studies of 
the Z boson mass peak and extrapolated to the signal re- 
gion. The total number of estimated background events 
is constrained to be the observed number. The di-photon 
invariant mass distribution for the events passing the full 
selection is shown in Fig. [2] The full- width at half maxi- 
mum of a signal with nih=120 GeV would be 4.2 GeV. 

The expected signal yield, summing all production pro- 
cesses, and estimated background composition for a total 
integrated luminosity of 38 pb^ 1 are summarised in Ta- 
ble [T] A total of 99 events passing all selection criteria 
are observed in data in the di-photon mass range from 
100 GeV to 150 GeV. The background in this region is 
modelled by fitting an exponential function to the data. 
The signal peak is modelled by a Gaussian core portion 
and a power-law low-end tail [53]. Tails in the signal res- 
olution are modelled by a wide Gaussian component of 
small amplitude. No significant excess of events over the 
continuous background is found for any Higgs boson mass. 
The systematic uncertainty on the total signal acceptance 
is ±15%, where the dominant contributions come from 
photon identification (±11%) and photon isolation effi- 
ciencies (±10%). 



5 Search for H -»■ WW 

The search for the Higgs boson in the decay channel H —> 
WW benefits from the large branching ratio of the Higgs 
boson to decay into a pair of W bosons for masses above 
ma ^ HO GeV, the sizable W boson decay rates to 
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Table 1. The number of expected and observed events in the 
H — > 77 search in the di-photon mass range from 100 GeV to 
150 GeV for an integrated luminosity of 38 pb _1 . Also shown 
is the composition of the background expected from Monte 
Carlo simulation and the division of the observed data as dis- 
cussed in the text as well as the expected number of H — > 77 
signal events for a Higgs boson mass of mn = 120 GeV. To- 
tal uncertainties are shown in the middle column while in the 
rightmost column the statistical and systematic uncertainties, 
respectively, are given. 





Expected 


Observed or Estimated 


Total 


120 ±27 


99 


77 


86±23 


75.0±13.3+36 


13 


31±15 


19.6±7.5±3.9 


33 


1±1 


1.5±0.7±os 


Z/l* 


2.7±0.2 


2.9±0.1±0.6 
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Fig. 2. Distribution of the di-photon invariant mass for the 
99 events from data passing all event selection criteria in the 
H — > 77 search and for the Monte Carlo prediction. The overall 
uncertainty on the expected total yield is illustrated by the 
yellow band. The uncertainty due to the reducible background 
is also shown (dark yellow band). The predictions for the main 
components of the background (di-photon, photon-jet, jet-jet 
and Drell-Yan) are also illustrated. 



leptons and the powerful identification of leptons with 
the ATLAS detector. It offers the greatest sensitivity of 
any search channel when the Higgs boson mass is close 
to twice the W boson mass, to# ~ 165 GeV. Two dif- 
ferent decay modes of the W bosons are considered: the 
H — > WW^-*' — > Ivlv channel is pursued for Higgs boson 
masses in the range from 120 GeV to 200 GeV, and the 
H — > WW — > Ivqq decay mode is used for Higgs boson 
masses in the range from 220 GeV to 600 GeV. The anal- 
yses are described below and further details can be found 
in Refs. [H1IS] 



5.1 Search for H -> WW^ -> Ivlv 

The H — > WW^ — > ivlv analysis is performed using 
a dataset corresponding to an integrated luminosity of 
35 pb -1 . Events are selected requiring exactly two iso- 
lated leptons with opposite charge. The leading lepton is 
required to have pt > 20 GeV and the sub-leading lep- 
ton is required to have pt > 15 GeV. Events are classified 
into three channels depending on the lepton flavours: e/j, 
ee or \x\x. If the two leptons are of the same flavour, their 
invariant mass (mu) is required to be above 15 GeV to 
suppress background from T production. To increase the 
sensitivity, the selections are then allowed to depend on 
the Higgs boson mass hypothesis. For all lepton combina- 
tions in the low (high) mass Higgs boson search, mu is 
required to be below 50 (65) GeV for Higgs boson masses 
rriH < 170 GeV (mH > 170 GeV) which suppress back- 
grounds from top-quark production and Z boson produc- 
tion. The missing transverse energy in the event is required 
to be E™ lss > 30 GeV. An upper bound is imposed on the 
azimuthal angle between the two leptons, A<fru< 1.3 (1.8) 
radians, taking advantage of the spin correlations |65) ex- 
pected in the Higgs boson decay. The signal region is de- 
fined by the transverse mass (ttit) |66j : 



TflT 



(£« 



£™ iss ) 2 



where E% 



|pmxss| 



P™ ss ) 2 , (1) 



#f 1ISS and 



is the transverse momentum of the dilepton system. The 
transverse mass is required to be 0.75-m# < tot < tuh for 
the event to be considered in a given mu range. Events are 
also treated separately depending on whether they have 
zero jets (0-jet channel) or one jet (1-jet channel) recon- 
structed with 1 77 1 < 4.5 due to the differences in back- 
ground composition and expected signal-to-background 
ratio. To suppress background from top-quark production, 
events in the 1-jet channel are rejected if the jet is iden- 
tified as coming from a 6-quark. Events with two or more 
jets have been analysed as a separate channel. However, 
due to the marginal contribution to the overall sensitiv- 
ity given the current total integrated luminosity and the 
additional systematic uncertainties, this channel is not in- 
cluded in this combination. 

The expected background contributions from WW, 
top-quark and VF+jets production are normalised using 
dedicated control regions in data as described in the next 
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sections. Other smaller backgrounds are normalised ac- 
cording to their theoretical cross sections. The background 
from Z/7*+jets production is normalised to the theoreti- 
cal cross section with a correction factor determined from 
data. 



5.1.1 The WW background 

The di-boson WW continuum can be distinguished from 
the Higgs boson signal through the kinematic selections. 
A control region is defined by changing the cut on mu to 
require over 80 GeV(but not within 10 GeV of the Z boson 
mass if the leptons are of the same flavour) and removing 
the selections on niT and A(f>u- The expected ratio of the 
background contribution in the control region and in the 
signal region is taken from Monte Carlo simulation. The 
three main sources of systematic uncertainty affecting this 
ratio are the theoretical uncertainty on the extrapolation, 
the jet energy scale uncertainty and the limited statistics 
in the simulated sample. Uncertainties due to these effects 
of ±6% in the 0-jet channel and ±17% in the 1-jet channel 
have been determined. 



5.1.2 The tt and single top-quark backgrounds 

Top-quarks, whether from strong interaction tt produc- 
tion or weak interaction single top-quark production, are 
a copious source of final states with one or two W bosons 
accompanied by one or more jets. Due to kinematic se- 
lection one or more of these jets may fail identification, 
thereby leading to a final state similar to that from the 
H WW signal. 

The background from top-quark production in the 0- 
jet channel is estimated by first removing the jet veto. 
This gives a sample dominated by top-quarks, and the ex- 
pected contamination from other processes in the control 
region is subtracted from the observed event yield. Then 
the probability that top events pass the jet veto is derived 
from the measured probability of not reconstructing a jet 
in data, using a sample of top candidates with two lep- 
tons, one b-jet and no other jet. The dominant systematic 
uncertainties originate from the limited statistics in data 
and the jet energy scale. A total uncertainty of ±60% has 
been determined for the top-quark background estimate 
in the 0-jet channel. 

The top-quark background in the 1-jet channel is nor- 
malised using a control region where the veto on jets com- 
ing from 6-quarks is reversed and the A<pu, mil and mx 
selections are removed. An extrapolation factor from the 
control region to the signal region is estimated from Monte 
Carlo simulation. The dominant systematic uncertainties 
on the top-quark background estimate in the 1-jet chan- 
nel are ±23% from the theoretical uncertainties on the 
extrapolation factor and ±22% from the uncertainty on 
the 6-tagging efficiency. 



5.1.3 The TV+jets background 

The production of W bosons accompanied by jets can 
mimic the H — > WW signal if one of the jets is mis- 
identified as an isolated lepton. The TV±jets background 
is normalised using a control region defined by relaxing the 
identification and isolation criteria for one of the two lep- 
tons. The contribution to the signal region is estimated 
by multiplying the rate measured in the control region 
by the probability for fake leptons which pass the relaxed 
identification and isolation criteria to also pass the original 
lepton selection criteria. This misidentification probability 
is measured in a multi-jet data sample. The major sources 
of systematic uncertainty for the W-^+jets background es- 
timate come from the bias introduced by the jet trigger 
threshold used to select the multi-jet events and the resid- 
ual difference in kinematics and flavour composition of the 
jets in multi-jet events and in events from M^+jets pro- 
duction. The total uncertainty on the estimated W / ±jets 
background is ±50%. 

5.1.4 The Z/7*+jets background 

The largest cross section for producing two isolated, high- 
Pt leptons comes from the Z/j* — > tt process. The back- 
ground from Z/7*+jets is significantly reduced by the up- 
per bound on mu and the requirement of high _E™ 1SS in 
the signal region. To correct for potential mis-modelling 
of the distribution of E™ lss at high values, a correction 
factor is derived from the observed difference between the 
fraction of events passing the E™ lss > 30 GeV selection 
in data and Monte Carlo simulation for events with m« 
within 10 GeV of the Z boson mass [H]. As the discrep- 
ancy between data and Monte Carlo tends to be larger in 
events with jets, the correction factor is larger in the 1-jet 
channel than in the 0-jet channel. The flavours of the two 
leptons in the event also impact the magnitude of the cor- 
rection factor, since any discrepancies between data and 
simulation have different sources. In the 1-jet channel, the 
correction factors are found to be 1.2±O.4±O.10in the ee 
analysis and 2.4±0.5±0.2 in the /i/i analysis. Under the as- 
sumption that the same correction factors apply to events 
below the upper bound on mu, the expected Z/7*+jets 
background is obtained from the Monte Carlo simulation 
normalised to the product of the theoretical cross section 
and the correction factors. 



5.1.5 Results for the H -s- WW^ Ivlv search 

The expected and observed numbers of events in the H — > 
WW analysis for a Higgs boson mass of 170 GeV are 
shown in Table [5J Three events in total are observed in 
the 0-jet channel for the combined ee, e/x and \x\x final 
states, compared to an expected number of events from 
background sources only of 1.70 ±0.12 ±0.17. More events 

3 When two errors are quoted the first is statistical and the 
second systematic. 
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are expected in the /i/i channel compared to the ee channel 
due to different lepton identification efficiencies for elec- 
trons and muons. In the 1-jet channel, one event is ob- 
served in the data compared to a total number of expected 
events from background sources of 1.26 ± 0.13 ± 0.23. The 
observed wit distributions in data after all selections ex- 
cept the transverse mass cut for the combined e/x, ee and 
\x\x channels are compared to the expected distributions 
from simulated events in Fig. [3] 
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Fig. 3. Distributions of the transverse mass ibt in the 0-jet 
channel (a) and 1-jet channel (b) for the H — > WW^ — > ivtv 
search after all selections except the transverse mass cut for 
the combined e/j,, ee and fifi channels. The error bars reflect 
Poisson asymmetric errors. A Higgs boson signal is shown for 
mn = 170 GeV. The selections applied for run = 170 GeV are 
indicated by the two vertical dotted lines. 



5.2 Search for H WW -> Ivqq 

The H — > WW — > Ivqq analysis uses a dataset corre- 
sponding to an integrated luminosity of 35 pb _1 . Events 



are selected requiring exactly one lepton with > 30 GeV. 
The missing transverse energy in the event is required to 
be E™ 1SS > 30 GeV. Events with fewer than two jets are 
rejectecQ- Events with > 4 jets are treated as a separate 
search channel, which is however not included in the cur- 
rent combination. The pair of jets with invariant mass 
closest to the W boson mass is considered to be coming 
from the W boson and the measured mass must be be- 
tween 71 GeV and 91 GeV. The event is rejected if any 
of the jets in the event is identified as coming from a b- 
quark. The invariant mass of the Higgs boson candidate, 
mivqqi is reconstructed with a W boson mass constraint 
on the lepton-neutrino system giving rise to a quadratic 
equation. If there are two solutions the one corresponding 
to the lower longitudinal momentum is taken; if complex 
the real part is used. 

The dominant source of background events in the H — > 
WW — > Ivqq search comes from VF+jets production. The 
contribution from QCD events is estimated by fitting the 
observed E™ lss distribution as the sum of templates taken 
from simulation. Table [3] shows the expected numbers of 
signal and background events in the signal region, as well 
as the observation. In the channel with only two and no 
additional jets, 450 events are observed in the data passing 
all selection criteria compared to an expected yield from 
background sources of 450±41 events. In the channel with 
one extra jet 263 events are observed, compared to an 
expected number of background events of 224 ± 15. 

The distributions of the invariant mass for the Higgs 
boson candidates in data are compared to the expected 
distributions from simulated events in Fig. 2) The mi uqq 
background spectrum is modelled with a falling exponen- 
tial function. The impact of the functional form has been 
investigated by replacing the single exponential with a 
double exponential, by histograms taken from simulation, 
and by a mixture of both methods without significant 
change in the results. It should be noted that the limit 
extraction is made using the exponential fit, not by com- 
parison with the simulated background. 

6 Search for H ZZ& 

Three different H — > ZZ^ final states are considered 
here: H ZZ^ till, H -> ZZ -> iivv and H -> 
ZZ -> liqq. In the H -> ZZ^ -> till search, the excel- 
lent energy and momentum resolutions of the ATLAS de- 
tector for electrons and muons lead to a narrow expected 
four-lepton invariant mass peak on top of a continuous 
background. The dominant background component is the 
irreducible ZZ& -> till process. In the low Higgs boson 
mass region, where one of the Z bosons is off-shell and de- 
cays into a pair of low transverse momentum leptons, the 
reducible backgrounds from Z+jets production and tt pro- 
duction are also important. For Higgs boson masses above 
rriH ^ 200 GeV both Z bosons are on-shell. In this region 
the decay modes H -> ZZ -> liqq and H — > ZZ — >• Uvv, 

4 In this channel the jet pt threshold is raised from 25 GeV 
to 30 GeV. 
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Table 2. Numbers of expected signal (rriH = 170 GeV) and background events and the observed numbers of events in the 
data passing all selections in the H — > WW^*' — > IvLv search. The dataset used in this analysis corresponds to an integrated 
luminosity of 35 pb . The uncertainties shown are the statistical and systematic uncertainties respectively. 



0-jet channel 





e/i 


ee 


/i/i 


WW 


0.71 ±0.05 ±0.06 


0.20 ±0.03 ±0.02 


0.53 ±0.02 ±0.05 


ti and single top 


0.09 ± 0.05 ± 0.06 


0.03 ±0.01 ±0.02 


0.08 ±0.04 ±0.06 


WZ/ZZ/W-y 


0.020 ±0.001 ±0.001 


0(< 0.001) ±0 


0.010 ±0.001 ±0.001 


Z/ 7 *±jets 


0(< 0.001) ±0 


0(< 0.001) ±0 


0(< 0.002) ±0 


W+jets 


0.01 ±0.01 ±0.01 


0.02 ±0.01 ±0.01 


0± 0.10 ±0.01 


Total Background 


0.83 ±0.07 ±0.13 


0.25 ±0.08 ±0.04 


0.62 ±0.05 ±0.10 


H -> WW W -> tulv 


0.62 ±0.01 ±0.18 


0.20 ±0.01 ±0.07 


0.44 ± 0.01 ±0.12 


Observed 


1 


1 


1 


1-jet channel 




e/i 


ee 


/i/i 


WW 


0.18 ±0.03 ±0.03 


0.05 ±0.02 ±0.01 


0.16 ±0.03 ±0.02 


ti and single top 


0.26 ±0.07 ±0.11 


0.10 ±0.02 ±0.04 


0.15 ±0.04 ±0.07 


WZ/ZZ/W-y 


0.01 ±0.001 ±0.001 


0(< 0.001) ±0 


0(< 0.001) ±0 


Z/ 7 *+jets 


0(< 0.01) ±0 


0.05 ±0.02 ±0.02 


0.25 ±0.08 ±0.05 


VK+jets 


0.02 ±0.02 ±0.01 


0.03 ±0.20 ±0.01 


0± 0.10 ±0.01 


Total Background 


0.47 ±0.08 ±0.16 


0.23 ±0.04 ±0.06 


0.56 ±0.09 ±0.14 


H -> WW { * } -> tulv 


0.31 ±0.01 ±0.09 


0.08 ±0.01 ±0.03 


0.21 ± 0.01 ±0.06 


Observed 








1 



Table 3. Numbers of expected signal (ran = 400 GeV) and background events and the observed numbers of events in the data 
passing all selections in the H — > WW — > ivqq search. The dataset used corresponds to an integrated luminosity of 35 pb -1 . 
The quoted uncertainties are combinations of the statistical and systematic uncertainties. 



ra H = 400 GeV 


H + 0-jets 


H + 1-jet 




euqq 


fivqq 


evqq 


[ivqq 


W/Z+jets 


157 ± 22 


259 ± 34 


39.1 ± 6.2 


119 ± 12 


Multi-jet 


11.1 ± 1.6 


4.5 ±0.6 


17.7 ± 2.8 


13.3 ± 1.3 


Top 


5.3 ± 1.7 


7.7 ±2.5 


15.5 ±5.0 


18.2 ±5.8 


Di-boson 


1.8 ±0.3 


3.0 ±0.4 


0.6 ±0.1 


0.9 ±0.1 


Total Background 


175 ± 22 


275 ± 34 


72.9 ±8.4 


151 ± 13 


H -> WW ivqq 


0.5 ±0.2 


0.6 ±0.2 


0.5 ±0.2 


0.5 ±0.2 


Observed 


177 


273 


87 


176 



which have substantially larger branching ratios but also 
larger backgrounds compared to the H — > ZZ^ — > till 
decay, provide additional sensitivity. The analyses of the 
H -> ZZ — S- ilqq and H -> ZZ — S- livv channels require 
that both Z bosons are on-shell, which limits the contribu- 
tion from the reducible backgrounds from Z±jets produc- 
tion and tt production. In this paper the H — > ZZ — > ilqq 
and H — > ZZ — > ilvv search channels have been used 
for Higgs boson masses in the range 200 GeV < m# < 
600 GeV, a range extending beyond the sensitivities of 
LEP and Tevatron experiments [HE]. Further details of 
the three analyses can be found in Refs. [T0llTT| . 

6.1 Search for H -> ZZ& -> lilt 

Candidate events are selected requiring two same-flavour 
and opposite-charge pairs of leptons. Muons with j»t > 



7 GeV and electrons with pt > 15 GeV are considered, 
while at least two out of the four leptons must satisfy 
Pt > 20 GeV. All leptons are required to be well sepa- 
rated from each other, isolated from other activity in the 
tracking detectors and the calorimeters and have low track 
impact parameters with respect to the primary vertex. At 
least one of the lepton pairs is required to have an invari- 
ant mass within 15 GeV (within 12 GeV if the combined 
four-lepton mass is high) of the Z boson mass. The re- 
quirement on the invariant mass of the second lepton pair 
varies as a function of the Higgs boson candidate mass, 
muu- The effective Higgs boson candidate mass resolu- 
tion a (ran), including the intrinsic width at the Higgs 
boson mass hypothesis being tested, is used to define an 
allowed range for the reconstructed Higgs boson candidate 
mass. The latter is required to be within ±5cr(m^f ) of the 
tested Higgs boson mass for the event to be considered. 
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(a) extra jet channel 
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(b) 1 extra jet channel 

Fig. 4. Distributions of the invariant mass mt vqq for the 
H — s> WW — > Ivqq search after the application of all selection 
criteria and the W-mass constrained fit. The background fit is 
shown as a continuous line. In (a) no extra jets are allowed and 
in (b) one additional jet is required. The Higgs boson signal is 
shown for mn = 400 GeV and the expected yield is scaled up 
by a factor of 30 for illustration purposes. 



The magnitude of the ZZ^ background is normalised 
to the measured Z boson cross section multiplied by the 
expected ratio of the cross sections <Jzz l&z from theoret- 
ical calculations [57] ■ This estimate is independent of the 
luminosity uncertainty, and the cross section ratio is less 
affected by theoretical uncertainties than the <izz cross 
section alone. The total uncertainty on the ZZ 1 -*^ back- 
ground estimate is ±15%. The reducible Z+jets back- 
ground arises predominantly from Z boson production in 
association with a pair of heavy flavour quarks which de- 
cay semi-leptonically. This background is normalised using 
dedicated control regions in data where the lepton iden- 
tification requirements are relaxed for the second pair of 
leptons. The final uncertainty on the Z+jets background is 
±20%. The ti background is estimated from Monte Carlo 
simulation and normalised to its theoretical cross section. 



A total uncertainty of ±25% is estimated for the ti back- 
ground contribution. 

After the application of all selection criteria, no can- 
didate events remain in data for the H -> ZZ^> -> Utt 
search at any Higgs boson mass. This is consistent with 
the small background and signal yields expected with the 
integrated luminosity of 40 pb -1 used in this analysis. 
The results are shown in Table 2] for two selected Higgs 
boson masses of ran = 130 GeV and mn = 200 GeV. The 
distribution of the Higgs boson candidate invariant mass, 
muu, before applying the lepton impact parameter and 
isolation requirements is shown in Fig. [SJ 



Table 4. Expected signal and background event yields in the 
H — > ZZ^ — > till search within ±5a(m,fr) for two selected 
Higgs boson masses. No events are observed in the data. The 
dataset used corresponds to a total integrated luminosity of 
40 pb -1 . The quoted uncertainties are combinations of the sta- 
tistical and systematic uncertainties. 



m H (GeV) 




130 


200 


Total backg 


ground 


0.010±0.002 


0.090±0.014 


H -> ZZ W 


-fim 


0.015±0.003 


0.095±0.017 


Observed 












100 150 200 250 300 350 400 450 500 550 600 

m NII [GeV] 



Fig. 5. Distribution of muu in the H -> ZZ^ -> ££££ search 
before applying the lepton impact parameter and isolation re- 
quirements which remove the two candidates. The error bars 
reflect Poisson asymmetric errors. 



6.2 Search for H -> ZZ -> Zlqq 

Events are selected requiring exactly two same-flavour lep- 
tons with an invariant mass 76 GeV < mu < 106 GeV 
and at least two jets. To reduce background from top pro- 
duction the missing transverse energy is required to be 
E™ iss < 50 GeV. The two jets in the event with the high- 
est individual j»t are required to have an invariant mass, 
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rrijj, in the range 70 GeV < nijj < 105 GeV. Additional 
background rejection is obtained for high mass by using 
the fact that the final state jets and leptons are boosted 
in the directions of the two Z bosons. For Higgs boson 
searches at win > 360 GeV, the two jets are required to 
have pt > 50 GeV. Furthermore, the azimuthal angles 
between the two jets, A<pjj, and between the two leptons, 
A(f>u, must both be less than tt/2. 

The Higgs boson candidate mass is constructed from 
the invariant mass of the two leptons and the two jets in 
the event, mujj- The two jets are constrained to have an 
invariant mass equal to the Z boson mass to improve the 
Higgs boson candidate mass resolution. 

6.2.1 Background estimates for the H — » ZZ — > Uqq 
search 

The dominant background in the H —> ZZ —> Uqq search 
channel is expected to come from Z+jets production. Other 
significant sources are ti production, multi-jet production 
and ZZ/WZ production. All backgrounds, except for the 
multi-jet background, are estimated from Monte Carlo 
simulation. For the Z+jets and the ti backgrounds the pre- 
dictions from simulation are compared against data in con- 
trol samples which are dominated by these backgrounds. 
The Z+jets control region is defined by modifying the 
rrijj selection to instead require 40 GeV < rrijj < 70 GeV 
or 105 GeV < rrijj < 150 GeV. The ti control region 
is denned by reversing the E!f lss selection and modifying 
the m« selection to require 60 GeV < m« < 76 GeV or 
106 GeV < m u < 150 GeV. Both the Z+jets and the 
tt background estimates from Monte Carlo simulation are 
found to be in good agreement with data in the control 
samples. The contribution from W+jets is very small and 
assumed to be adequately modelled. The multi-jet back- 
ground in the electron channel is derived from a sample 
where the electron identification requirements are relaxed. 
In the muon channel, the multi-jet background is taken 
from Monte Carlo after verifying the accuracy of the sim- 
ulation using a data sample where the two muons in the 
event are required to have the same charge. 

6.2.2 Results for the H — s- ZZ -> Uqq search 

The H — > ZZ — > Uqq analysis is performed for Higgs 
boson masses between 200 GeV and 600 GeV in steps of 
20 GeV. Table [5] summarises the numbers of estimated 
background events and observed events in data for the se- 
lections below and above mn — 360 GeV. The numbers of 
expected signal events for two representative Higgs boson 
masses are also shown. For the low mass search, 216 events 
are observed in data passing all selection criteria compared 
to an expected number of events from background sources 
only of 226±4±28 events. The corresponding numbers for 
the high mass searches are 11 events observed in data com- 
pared to an expected yield of 9.9 ± 0.9 ± 1.5 events from 
background sources only. The distribution of the recon- 
structed Higgs boson candidate mass mujj for the events 
passing all of the selection criteria is shown in Fig. [51 




100 200 300 400 500 600 700 800 900 1000 



miijj [GeV] 

Fig. 6. Distribution of miijj for events passing all of the se- 
lection criteria in the H — > ZZ — > Uqq search. The expected 
yield for a Higgs boson with a mass ran = 300 GeV is also 
shown, multiplied by a factor of 20 for illustrative purposes. 
The contribution labelled "Other" is mostly from top events 
but includes also QCD multijet production. 



6.3 Search for H -> ZZ -> Uvv 

The H — > ZZ — > Uvv final state is characterised by two 
charged leptons and large E^ iss . Events are selected by 
requiring exactly two leptons of the same flavour with an 
invariant mass 76 GeV < m« < 106 GeV. Events are re- 
jected if any jet is identified as coming from a &-quark. 
The selection has been optimised separately for searches 
at low (m„ < 280 GeV) and high (m H > 280 GeV) val- 
ues of the Higgs boson mass. Events are required to have 
E™ iss > 66 (82) GeV and A<f» u < 2.64 (2.25) radians 
for the low (high) mass region. For the low mass region 
A(j)u > 1 radian is also required. The Higgs boson candi- 
date transverse mass is obtained from the invariant mass 
of the two leptons and the missing transverse energy. 

6.3.1 Background estimates for the H — s- ZZ — > Uvv 
search 

A major background in the H — > ZZ — ¥ Uvv search chan- 
nel comes from di-boson production and is estimated from 
Monte Carlo simulation. Background contributions from 
tt and W^+jets production are also obtained from Monte 
Carlo simulation, and the estimated yields arc verified by 
comparing with the number of observed events in dedi- 
cated control samples in the data. Both the tt and the 
VF+jets control regions are defined by modifying the mu 
selection to instead require 60 GeV < mu < 76 GeV or 
106 GeV < mu < 150 GeV. The tt control region also 
requires that the events pass E!f ms > 20 GeV and that at 
least one jet is identified as coming from a 6-quark. The 
W^+jets control region instead requires _E™ ISS > 36 GeV 
and that no jets in the events are identified as coming 
from a 6-quark. The observed event yields in the con- 
trol regions for ti and VF+jets production are in good 
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Table 5. Numbers of events estimated as background, observed in data and expected from signal in the H — > ZZ — s> llqq search 
for low mass {tuh < 360 GeV) and high mass (run > 360 GeV) selections. The signal, quoted at two mass points, includes small 
contributions from lilt and llvv decays. Electron and muon channels are combined. The uncertainties shown are the statistical 
and systematic uncertainties, respectively. 



Source 


low mass selection 


high mass selection 


Z+jets 


214±4±27 


9.1 ±0.9 ± 1.4 


VK+jets 


0.33 ±0.16 ±0.17 




ti 


0.94 ±0.09 ±0.25 


0.08 ±0.02 ±0.03 


Multi-jet 


3.81 ±0.65 ±1.91 


0.11 ±0.11 ±0.06 


ZZ 


3.80 ±0.10 ±0.73 


0.30 ±0.03 ±0.06 


WZ 


2.83 ±0.05 ±0.88 


0.29 ±0.02 ±0.10 


Total background 


226 ± 4 ± 28 


9.9 ± 0.9 ± 1.5 


H -> ZZ llqq 


0.60 ± 0.01 ± 0.12 (m H = 200 GeV) 


0.24 ± (< 0.001) ± 0.05 (m H = 400 GeV) 


Observed 


216 


11 



agreement with the predictions from the Monte Carlo sim- 
ulation. The background from Z±jets production is es- 
timated from Monte Carlo simulation after comparison 
studies of the £"™ ISS distribution between Monte Carlo and 
data. The multi-jet background in the electron channel is 
derived from a sample where the electron identification 
requirements are relaxed. In the muon channel, the multi- 
jet background is estimated from a simulated sample of 
semi-leptonically decaying b- and c-quarks and found to be 
negligible after the application of the mu selection. This 
was verified in data using leptons with identical charges. 



6.3.2 Results for the H — s- ZZ — > llvv search 



The H ZZ llvv analysis is performed for Higgs 
boson masses between 200 GeV and 600 GeV in steps of 
20 GeV. Table [5] summarises the numbers of events ob- 
served in the data, the estimated numbers of background 
events and the expected numbers of signal events for two 
selected to# values. For the low mass selections, five events 
are observed in data compared to an expected number of 
events from background sources only of 5.8±0.5± 1.3. The 
corresponding results for the high mass selections are five 
events observed in data compared to an expected yield of 
3.5±0.4±0.8 events from background sources only. In ad- 
dition to the H — > ZZ — > llvv decays, several other Higgs 
boson channels give a non-negligible contribution to the 
total expected signal yield. In particular, H — > WW^ — > 
Ivlv decays can lead to final states that are very similar 
to H — > ZZ — > llvv decays. They are found to contribute 
significantly to the signal yield at low m# values. The 
expected number of events from H — > VFWW — > Ivlv de- 
cays relative to that from H — > ZZ — > llvv decays is 76% 
for m H = 200 GeV and 9% for m H = 300 GeV. The kine- 
matic selections prevent individual candidates from being 
accepted by both searches. The E™ lss distribution before 
vetoing events with low _B™ 1SS is shown in Fig. [7] 




E™ ss [GeV] 

Fig. 7. Distribution of missing transverse energy in the H — > 
ZZ — > llvv search in the electron channel before vetoing events 
with low _E™ 1SS . The expected yield for a Higgs boson with 
uih = 400 GeV is also shown. The distribution in the muon 
channel is similar with four events seen which have _E™ 1SS above 
80 GeV. 



7 Combination method 

The limit-setting procedure uses the power-constrained 
profile likelihood method known as the Power Constrained 
Limit, PCL [T21[TU[SB]. This method is preferred to the 
more familiar CL S |15j technique because the constraint 
is more transparently defined and it has reduced overcov- 
erage resulting in a more precise meaning of the quoted 
confidence level. The resulting PCL median limits have 
been found to be around 20% tighter than those obtained 
with the CL S method in several Higgs searches. The ap- 
plication of the PCL method to each of the individual 
Higgs boson search channels is described in Refs. [TMll|. 
A similar procedure is used here. The individual analyses 
are combined by maximising the product of the likelihood 
functions for each channel and computing a likelihood ra- 
tio. A single signal normalisation parameter fi is used for 
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Table 6. Numbers of events estimated from background, observed in data and expected from signal in the H — > ZZ — > llvv 
search for low mass (run < 280 GeV) and high mass (rriH > 280 GeV) selections. Electron and muon channels are combined. 
The expected signal events include minor additional contributions from H — > ZZ — > iiqq , H — > ZZ^ — > iiii and one which 
can be large from H — > WW^*' — > ivlv. The uncertainties shown are the statistical and systematic uncertainties, respectively. 



Source 


low mass selection 


high mass selection 


Z+jets 


1.09 ±0.29 ±0.59 


1.01 ±0.29 ±0.58 


W+jets 


1.07 ±0.31 ±0.64 


0.41 ±0.19 ±0.22 


tt 


1.90 ±0.10 ±0.63 


0.91 ±0.07 ±0.31 


Multi-jet 


0.11 ±0.11 ±0.06 




ZZ 


0.58 ±0.01 ±0.11 


0.51 ±0.01 ±0.10 


WZ 


0.57 ±0.01 ±0.10 


0.45 ±0.01 ±0.09 


WW 


0.43 ±0.02 ±0.09 


0.16 ±0.01 ±0.04 


Total background 


5.8 ± 0.5 ± 1.3 


3.5 ±0.4 ±0.8 


H -> ZZ -> llvv 


0.19 ± (< 0.001) ± 0.04 (m H = 200 GeV) 


0.30 ± (< 0.001) ± 0.06 (toh = 400 GeV) 


Observed 


5 


5 



Table 7. Summary of systematic uncertainties (in percent) of the signal yield. The correlated systematic uncertainties are 
given in detail, the uncorrelated ones are lumped together. The uncertainties are evaluated for a Higgs boson mass of 115 GeV 
for the ff^77 channel, 160 GeV for H -> WW M -> ivlv, 200 GeV for the H -> ZZ M -> till and 400 GeV for the 
remaining channels. Systematic errors marked with a dash are neglected. In the three channels, the impact of the lepton energy 
scale and resolution uncertainties on the efficiency was found to be negligible, but they can still influence the fit via the signal 
distributions. In the H — > WW —¥ ivqq channel a jet energy scale uncertainty can only decrease the efficiency; the resolution 
uncertainty is negligible in comparison. 





77 


H -> 


WW 




H -+ ZZ 








ivlv 


Ivqq 


llll 


iivv 


llqq 


Luminosity 


±3.4 


±3.4 


±3.4 


±3.4 


±3.4 


±3.4 


e/7 efficiency 


± 11 


±8.2 


±2.6 


±2.1 


±4.6 


±0.0 


e/7 energy scale 




±1.6 






±0.5 


±0.2 


e/7 resolution 




±1.6 






±0.1 


±0.1 


/i efficiency 




±0.5 


±1.0 


±0.8 


±0.0 


±2.0 


fj, energy scale 




±4.8 






±1.2 


+0.2 
-2.2 


fj, resolution 




±1.2 






±0.1 


±0.1 


Jet energy scale 




±3.7 


-26 




±0.4 


+ 2.9 

-7.0 


Jet energy resolution 










±0.2 


+0.0 
-1.3 


6-tag efficiency 










±0.4 




Uncorrelated 


±10 


±5.0 











all analyses, where fj, is the ratio of the hypothesised cross 
section to the expected Standard Model cross section. 

Each channel has sources of systematic uncertainty, 
some of which are common with other channels. Table [7] 
lists the common sources of systematic uncertainties, which 
are taken to be 100% correlated with other channels. Let 
the search channels be labelled by I (I = H — >• 77, H — s- 
WW, ■ ■ ■), the background contribution, j, to channel I 
by ji and the systematic uncertainties by i (i = luminos- 
ity, jet energy scale, . . . ). The relative magnitude of the 
effect on the Higgs boson signal yield in channel I due to 
systematic uncertainty i is then denoted by ef i; and on 
background contribution ji, e^ h . The e/i's are constants; 
an individual can be zero if the channel in question 
is not affected by this source of systematic uncertainty. 
A common systematic uncertainty, i, which is shared be- 
tween channels I and V implies that and ef, i are both 
different from zero. If a systematic source i is shared be- 



tween the signal in channel I and background contribution 
ji then both e a u and e^ h are non-zero. For each source of 
systematic uncertainty i there is a corresponding nuisance 
parameter Si and an associated auxiliary measurement mj 
on a control sample (e.g. sidebands in a mass spectrum) 
that is used to constrain the parameter. The Si and mi 
are scaled so that Si = corresponds to the nominal ex- 
pectation and Si = ±1 corresponds to the ±ler variations 
of the source. When constructing ensembles for statistical 
evaluation, each rrii is sampled according to G(mi\Si, 1), 
the standard normal distribution. Using this notation, the 
total number of expected events in the signal region for 
channel I is given by: 

Nf xp = nL ui + Aih) + 

i 

E^IK 1+ 4^) ( 2 ) 
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for luminosity L, Standard Model cross sections ai (in- 
cluding efficiencies and acceptances), and expected back- 
grounds bji. Background estimates bji may come either 
from Monte Carlo simulations or from control regions in 
which the expected number of events, hji, is proportional 
to the expected background, via bji = ajifiji. Given the 
number of observed events in the signal region N° bs , the 
likelihood function can be written as: 



Pois{N? bs \Nf xp ) 



JJPois(n3-,|n i i)II G! ( m *l*' 1 ) 



where riji are the observed numbers of background events 
in the control regions and Pois(a;|y) is the Poisson proba- 
bility of observing x events given an expectation y. 

The combined likelihood is given by the product of the 
individual likelihoods for each channel: 



where I is implicitly an index over the individual histogram 
bins within the channels that used a binned distribution 
of a discriminating variable. 
The profile likelihood ratio 



A( M ) 



£(m,9Gu)) 

£(M>)) 

£(0,9(0)) ' 



A>0, 

A < o, 



(3) 



is computed by maximising the likelihood function twice: 
in the numerator /x, the ratio of the hypothesised cross 
section to the expected Standard Model cross section, is 
restricted to a particular value and in the denominator \i 
is allowed to float. The set of all nuisance parameters Si 
and nji is denoted 9. The maximum likelihood estimates 

of /I and 8 are denoted fi and 8, while 6(fi) denotes the 
conditional maximum likelihood estimate of all nuisance 
parameters with /i fixed. In this analysis the range of /i 
is restricted to the physically meaningful regime, i.e. it is 
not allowed to be negative. The test statistic is defined 
to be 



% = 



-21nA(/i), fx < fi, 
0, fi > fi, 



-2 In £ ("- 9 >» , fi < 0, 

£(0,(3(0)) 
C(p,,6 

0, fi > fi. 



< fi < M, 



(4) 



Monte Carlo pseudo-experiments are generated to con- 
struct the probability density function _/* ( g M | /j. , 8{p)) under 
an assumed signal strength /i, giving a p- value 



/(g /1 |/i,6'(/i))dq (U 



(5) 



To find the upper limit on /i at 95% confidence level, 

fi is varied to find Py. = 5%. Similarly, background-only 



Monte Carlo pseudo-experiments are used to find the me- 
dian fj, mec i along with the ±lcr and +2a bands expected 
in the absence of a signal. The procedure so far can be 
referred to as a CL s b limit. To protect against excluding 
the (signal) null hypothesis in cases of downward fluctu- 
ations of the background, the observed limit is not al- 
lowed to fluctuate below the —la expected limit. This is 
equivalent to restricting the interval to cases in which the 
statistical power of the test of /i against the alternative 
fi = is at least 16%. This is referred to as a Power Con- 
strained Limit. If the observed limit fluctuates below the 
16% power, the quoted limit is fi m ed — If. 



8 Systematic uncertainties in the combination 

The systematic uncertainty related to the luminosity is 
±3.4% and is fully correlated among all channels. It af- 
fects background estimates that are normalised to their 
theoretical cross sections; for most channels this is only 
true for backgrounds that are known to be small. In the 
H — s- ZZ -> ilvv and H -> ZZ -> Uqq channels ma- 
jor backgrounds are normalised to their theoretical cross 
sections, but in the latter case this is only done after com- 
paring with control regions. 

Sources of systematic uncertainty related to the event 
reconstruction are correlated between all the Higgs boson 
search channels. The uncertainty on the efficiency to re- 
construct electrons varies between 2.5% (central high-px 
electrons) and 16% (pt near 15 GeV, the lowest value used 
here) but it is assumed to be completely correlated. For 
muons the efficiency uncertainty ranges between 0.4% and 
2%. The jet systematic errors are typically larger for the 
channels where jets are explicitly required. They are dom- 
inated by the jet energy scale as the resolution effects tend 
to partially cancel and the i?™ ISS uncertainties are largely 
by-products of the uncertainties already discussed. 

The effect on the signal yield in each channel of the 
major sources of systematic uncertainty is summarised in 
Table [7J Uncertainties are treated as either uncorrelated 
or 100% correlated among channels. The largest uncorre- 
lated errors are photon isolation in H — ¥ 77 and jet rates 
in H — > WW(*> — > lvlv\ the latter is in principle cor- 
related with the H — > WW — > Ivqq channel but these 
channels are never used in the same mass region. Most 
backgrounds have been estimated by means of indepen- 
dent control samples; these estimates are assumed to be 
uncorrelated between the channels. 

Systematic uncertainties on the signal shape are ac- 
counted for in the H — > 77, H — > ZZ — > Uqq and 
H — > ZZ — > ilvv channels by considering three possible 
distributions and interpolating between them. Small sig- 
nal shape systematic uncertainties in the H — > ZZ^ — > 
(Mi and H — > WW — > ivqq channels are neglected. For 
m,H > 200 GeV the correlations in the shape systematics 
are taken into account and are treated as correlated with 
the signal normalisation uncertainties. 

The width of the Higgs boson signal at high mass is 
taken from the PYTHIA Monte Carlo [S3]. This underes- 
timates the width and the accepted cross section is con- 
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servatively scaled down by the ratio of the widths given 
in Ref. (69], which reached a maximum of 8% at 600 GcV, 
in all plots showing a ratio to the Standard Model. 

The systematic uncertainty coming from the total the- 
oretical Higgs boson cross section is not included in the 
combination and is shown separately in the figures as an 
uncertainty on the predicted cross section. 



9 Combination 

Each Higgs boson search channel is only sensitive for a 
range of Higgs boson masses. The ranges in which the var- 
ious channels have been analysed are detailed in Table 
In the H -> 77, H -> ZZ -> tlvv and H -> ZZ -> llqq 
channels, the final result is extracted from a fit of signal 
plus background contributions to the observed Higgs bo- 
son candidate mass distributions. In all other channels, 
limits are extracted from a comparison of the numbers 
of observed events in one or more signal regions to the 
numbers of estimated background events. 

Table 8. The Higgs boson mass regions in which individual 
search channels have been analysed. 





Mode 


Mass ranj 


ge, ( 




H — > 77 


110 - 


140 


H -> 


WW (t) -> ivlv 


120 - 


200 


H - 


-> WW -i ivqq 


220 - 


600 


H - 


-> ZZ^ -> mi 


120 - 


600 


H 


-> ZZ — > tlvv 


200 - 


600 


H 


-» ZZ — >• llqq 


200 - 


600 



The individual channels are shown in Fig. [3] in terms 
of the observed and the expected upper limits on (t/ctsm 
at the 95% confidence level. The step with which the limit 
is extracted, 5-10 GeV, does not match the H — > 77 reso- 
lution which has a full-width at half maximum of 4.4 GeV. 
However, it was established in Ref. [7| that no important 
fluctuations are missed. 

The search channels are grouped by the primary Higgs 
boson decay mode searched for, 77, WW or ZZ, and the 
limit on each mode is extracted in terms of the cross sec- 
tion for the process intended. Some channels have a con- 
tribution from signal modes other then the intended one. 
This is only significant for the H — » ZZ — > tlvv search, 
as discussed in Section [6.31 and implies that the ZZ limit 
assumes the Standard Model ratio between H to ZZ and 
H to WW decays. In addition, the WW search requires 
zero or one jet and is essentially designed for a spin-zero 
object produced largely via gluon fusion. The upper lim- 
its at 95% confidence level observed and expected in the 
absence of a signal are compared with the cross section 
expected for a Standard Model Higgs boson in Fig. |H1 

The combination of all channels is tested and the v^=Q 
in these fits varies between 7% and 89%, which does not 
suggest the presence of a signal. The combination of all 



channels is shown in Fig. [TU] in terms of the observed and 
the expected upper limit at the 95% confidence level. The 
statistical accuracy of the toy Monte Carlo used to extract 
the limits is about 5% on the observed limits and 7% on 
the expected ones, with somewhat larger variation on the 
edges of the one and two a bands. 

The excluded signal strength as a function of mn is 
summarised in Table |H1 using the PCL method. The re- 
sults are also calculated using the CL S method for compar- 
ison purposes: the extracted limits from both procedures 
are shown in Fig. [TTJ Also given is the profile likelihood 
ratio of a Standard Model Higgs boson to background only 
and the consistency of the data with the background-only 
hypothesis, p M= o- 

The results have been interpreted in terms of the heavy 
mass fourth generation model introduced in Section [3] 
This involves rescaling the gluon fusion component of the 
production cross section and the Higgs boson decay branch- 
ing ratios. The limits are shown in Fig. 1121 This model is 
excluded for Higgs boson masses between 140 GeV and 
185 GeV, while the region in which exclusion might be 
expected is between 136 GeV and 208 GeV. 

The excluded cross section ratios are summarised in 
Table M 



10 Conclusions 

The ATLAS search for the Standard Model Higgs boson 
in the mass range from 110 GeV to 600 GeV in 35 to 40 
pb _1 of data recorded in 2010 has been presented. With 
this luminosity there is not sufficient sensitivity to exclude 
the Standard Model Higgs boson, nor is there any evidence 
of an excess of events over the predicted background rates. 
However, these results give the most stringent constraints 
to date for Higgs boson masses above 250 GeV and are 
close to the Tevatron limits [5] at intermediate masses. 

An extension of the Standard Model, assuming the 
Higgs mechanism and adding a fourth generation of heavy 
quarks and leptons, is excluded for Higgs boson masses 
between 140 and 185 GeV. 
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Fig. 8. The expected and observed cross section limits, normalized to the Standard Model cross section, as a function of the 
Higgs boson mass for the individual search channels. The visually most apparent difference between expected and observed is 
in the H — > WW — > luqq channel, which has a deficit approaching one sigma both at 320 GeV and 480 GeV. These results use 
the profile likelihood method with a power constraint (PCL). The limits are calculated at the masses marked with symbols. 
The lines between the points are to guide the eye. The grey horizontal bands show the uncertainty on the Standard Model cross 
section prediction, with the inner region highlighting the contribution of QCD scale uncertainties. 




m H [GeV] 

Fig. 9. The expected and observed 95% PCL limits on the total cross section of a particle produced like the Standard Model 
Higgs boson and decaying with the width predicted by PYTH1A [53] to pairs of bosons: 77, WW or ZZ . The limits are 
calculated at the masses marked with symbols. The lines between the points are to guide the eye. The coloured bands show the 
cross section predictions and their uncertainties, with the inner region highlighting the contribution of QCD scale uncertainties. 
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100 200 300 400 500 600 

m H [GeV] 

Fig. 10. The expected and observed upper limits on the total cross section divided by the expected Standard Model Higgs 
boson cross section. This is a 95% PCL limit. The green and yellow bands indicate the range in which the limit is expected 
to lie in the absence of a signal. The fine dotted line marks the results obtained using CL S (,, and the application of the power 
constraint gives the solid line. The limits are calculated at the masses marked with symbols and the lines between the points 
are to guide the eye. 
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Fig. 11. Same as Fig. 1101 except that limits calculated using the CL S procedure are added. As expected, when the observed 
limits fluctuate up, both methods converge, but downward fluctuations are less pronounced with CL S due to its larger over- 
coverage. The fine dotted line marks the results obtained using CL S (,, and the application of the power constraint gives the 
solid line. The limits are calculated at the masses marked with symbols. The lines between the points are to guide the eye. The 
regions excluded by the combined LEP experiments [4] and the Tevatron experiments [6] are indicated. 
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Table 9. The signal cross sections, in multiples of the Standard Model cross section, that are excluded, and expected to be 
excluded, at 95% CL. The expected variation at ±ler is also given for the PCL limits. The bold numbers show the limit which 
should be used; for mass of 500 and 520 GeV the power constraint is applied. The likelihood ratio of signal plus background 
to background is also shown, as is the p-value (modified to go between and 1) for /i = 0, which can be used to estimate the 
discovery significance. 



raH(GeV) PCL limits CL S limits -2 In ^2 p-values 

C(0,B) 





Obs. 


-la 


Median 


+la 


Obs. 


Median 




jDm— n 

ff-L — U 


110 


18.7 


6.9 


21.5 


38.2 


28.1 


29.6 


0.1 


0.58 


115 


42.4 


7.8 


20.9 


34.7 


43.5 


25.3 


-0.3 


0.07 


120 


18.2 


4.3 


11.4 


19.9 


19.7 


15.4 


-0.3 


0.22 


130 


10.0 


2.5 


6.4 


10.9 


11.0 


8.5 


-0.6 


0.23 


140 


5.0 


1.7 


4.3 


7.6 


6.1 


5.9 


0.0 


0.41 


150 


2.5 


1.3 


3.2 


5.5 


4.0 


4.4 


1.0 


0.65 


160 


1.6 


0.7 


2.3 


3.9 


2.8 


3.1 


1.6 


0.68 


170 


3.4 


0.8 


2.4 


3.9 


3.8 


3.1 


-0.4 


0.27 


180 


5.3 


1.0 


3.0 


5.4 


5.6 


4.2 


-0.9 


0.18 


190 


8.8 


1.9 


4.8 


7.8 


9.2 


6.3 


-1.1 


0.11 


200 


9.7 


2.1 


5.4 


9.5 


9.9 


7.5 


-1.1 


0.15 


220 


15.9 


2.9 


10.0 


15.4 


17.1 


12.9 


0.0 


0.13 


240 


8.3 


2.5 


9.1 


14.5 


11.2 


11.9 


0.3 


0.57 


260 


7.4 


2.9 


7.7 


12.6 


10.8 


10.9 


0.4 


0.56 


280 


10.0 


2.5 


7.3 


13.9 


11.5 


10.2 


0.2 


0.31 


300 


9.5 


1.9 


7.2 


13.2 


11.4 


10.1 


0.2 


0.32 


320 


7.1 


2.5 


6.2 


12.1 


9.8 


9.5 


0.4 


0.40 


340 


9.0 


2.8 


6.5 


11.9 


9.9 


9.6 


0.4 


0.28 


360 


5.5 


2.5 


7.7 


12.5 


8.5 


9.5 


0.4 


0.63 


380 


5.6 


2.2 


7.3 


12.9 


8.6 


9.5 


0.4 


0.53 


400 


7.5 


2.4 


7.8 


13.6 


9.4 


9.6 


0.1 


0.49 


420 


8.0 


2.9 


8.4 


14.7 


10.4 


10.4 


0.2 


0.46 


440 


9.8 


3.5 


9.9 


15.8 


11.7 


11.8 


0.1 


0.46 


460 


6.3 


3.1 


8.7 


18.1 


10.1 


12.3 


0.4 


0.64 


480 


5.6 


3.7 


10.1 


17.7 


10.4 


13.4 


0.4 


0.80 


500 


3.1 


5.0 


9.9 


19.3 


11.8 


15.9 


0.5 


0.89 


520 


4.8 


5.3 


14.1 


23.1 


13.9 


18.5 


0.4 


0.86 


540 


6.3 


5.6 


16.7 


26.0 


16.8 


21.3 


0.4 


0.82 


560 


8.0 


6.6 


16.9 


33.5 


19.3 


23.7 


0.3 


0.80 


580 


19.7 


10.2 


22.2 


37.9 


27.5 


28.8 


0.1 


0.56 


600 


26.1 


10.6 


24.0 


49.2 


34.4 


33.9 


0.1 


0.45 
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Fig. 12. Same as Fig. [TO] but comparing the excluded cross section to the expected one when a fourth generation of high mass 
quarks and leptons with Standard Model-like couplings to the Higgs boson are included in the cross section calculations. The 
arrows indicate the regions excluded by the CMS experiment [12] and the Tevatron experiments [70] . The point set at 500 GeV 
is at the limit allowed by the power constraint. The limits are calculated at the masses marked with symbols. The lines between 
the points are to guide the eye. 



Table 10. The signal cross sections, in multiples of the high mass fourth generation model [50] cross section, that are excluded, 
and expected to be excluded, in the absence of signal, at 95% CL. 



mjf(GeV) 




PCL limits 


CL S limits 


m H (GeV) 




PCL limits 


CL S limits 




Obs. 


-la 


Median 


Median 


Obs. 




Obs. 


-la 


Median 


Median 


Obs. 


110 


17.1 


7.8 


19.6 


25.2 


23.9 


320 


1.3 


0.4 


1.1 


1.7 


1.8 


115 


35.8 


6.7 


17.0 


22.4 


36.6 


340 


1.7 


0.5 


1.2 


1.9 


1.9 


120 


4.3 


1.2 


2.8 


3.8 


4.8 


360 


1.3 


0.6 


1.7 


2.0 


1.9 


130 


2.0 


0.5 


1.2 


1.6 


2.1 


380 


1.5 


0.7 


1.8 


2.1 


2.0 


140 


0.9 


0.2 


0.7 


0.9 


1.0 


400 


2.0 


0.8 


1.9 


2.2 


2.3 


150 


0.4 


0.2 


0.4 


0.6 


0.5 


420 


2.1 


0.9 


2.0 


2.6 


2.5 


160 


0.2 


0.1 


0.3 


0.4 


0.3 


440 


2.6 


1.1 


2.2 


2.8 


3.0 


170 


0.4 


0.1 


0.3 


0.4 


0.4 


460 


1.6 


1.2 


2.4 


3.0 


2.6 


180 


0.6 


0.1 


0.4 


0.5 


0.7 


480 


1.5 


1.0 


2.5 


3.3 


2.9 


190 


1.1 


0.2 


0.6 


0.8 


1.1 


500 


0.8 


1.4 


3.0 


3.5 


3.0 


200 


1.2 


0.3 


0.7 


1.0 


1.3 


520 


1.3 


0.9 


3.0 


4.1 


3.6 


220 


2.3 


0.5 


1.2 


1.8 


2.4 


540 


1.7 


1.5 


3.3 


4.3 


3.8 


240 


1.2 


0.4 


1.1 


1.7 


1.7 


560 


1.8 


1.8 


4.3 


5.4 


4.6 


260 


1.1 


0.4 


1.1 


1.6 


1.6 


580 


4.8 


1.8 


4.6 


6.0 


6.3 


280 


1.6 


0.4 


1.1 


1.6 


1.9 


600 


6.4 


2.5 


5.5 


7.4 


7.8 


300 


1.5 


0.4 


1.1 


1.6 


1.9 
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